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^J ' We consider a model of quantum computation in which the set of operations is limited to 

nearest-neighbor interactions on a 2D lattice. We model movement of qubits with noisy 
(__J ■ SWAP operations. For this architecture we design a fault-tolerant coding scheme using the 

>D ' concatenated [[7, 1, 3]] Steane code. Our scheme is potentially applicable to ion-trap and 

C_^ ' solid-state quantum technologies. We calculate a lower bound on the noise threshold for 

our local model using a detailed failure probability analysis. We obtain a threshold of 
1.85 X 10 -5 for the local setting, where memory error rates are one-tenth of the failure 
rates of gates, measurement, and preparation steps. For the analogous nonlocal setting, 
we obtain a noise threshold of 3.61 X 10 — 5 . Our results thus show that the additional 
SWAP operations required to move qubits in the local model affect the noise threshold 
only moderately. 
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1 Introduction 

In order to perform large-scale quantum computation reliably, any quantum computer tech- 
nology will require error correction. Error-correcting quantum data is useful if the logical 
noise level is lowered by coding and error-correcting. Since error correction is noisy, a de- 
crease in logical noise level only occurs for sufficiently low initial noise rates, below a certain 
noise threshold [H HJ OH H] . It is clear that it is of great importance for the future of quan- 
tum computation to determine the values of such noise thresholds since they set the level of 
accuracy to be achieved by experimental programs. 

Over the past 10 years estimates for noise thresholds for various code architectures [U El 
[3 [8] have been determined. These estimates, depicted in the overview in Fig. 1 of [9], have 



"kmsvoreOcs. columbia.edu 
kdivinccQwatson. ibm.com 
c bterhal@gmail.com 
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varying degrees of accuracy and rigor. In Rcf. [10j . entirely rigorous criteria for fault-tolerant 
design and threshold analysis were formulated. These rigorous criteria allow, at least for small 
error-correcting codes, an accurate determination of the noise threshold for simple stochastic 
noise models. 

In many threshold analyses it has been assumed that interactions between qubits can 
be performed in a manner independent of their spatial separation. In reality, interactions 
between stationary qubits are typically restricted to occur between nearest neighbors in some 
low-dimensional spatial layout, the spatial architecture. A transport mechanism will thus be 
required to bring qubits together. An example is the scalable ion-trap architecture where ions 
are shuttled between interconnected ion traps [11] . Transport mechanisms will also be required 
in solid-state implementations, such as the method involving electron spins in quantum dots 
[12] . phosphorus spin qubits in silicon [13j . Josephson-junction qubits [14] or neutral atoms 
in a 2D optical lattice [15] . 

Previous work [16j HZ] nas established the existence of a noise threshold in a local 
model. In Ref. [8] we made a first attempt to also estimate the effect of locality on the noise 
threshold, in particular the dependence on a code-dependent scale parameter and the failure 
rates of qubit transportation. In Rcf. [18], the authors derive a noise threshold of O(10 -7 ) 
for the [[7,1,3]] code for a quasi-onc-dimensional architecture. Some thresholds related to ion- 
trap systems are also analyzed in Refs. [HI [20]. For spin-qubits, Rcfs. [21] and [22] estimate 
noise thresholds taking into account transport mechanisms in conjunction with explicit two- 
dimensional layouts. None of these papers employ the rigorous methods of analysis outlined 
in Ref. [10] . We believe the use of rigorous methods can considerably impact the threshold 
values that are found. We briefly argue in Section [6J why the methodology laid out in Rcf. 
[10] gives fairly tight bounds for the fault-tolerance threshold. This implies that some of the 
previous threshold estimates for the [[7,1,3]] code have been too optimistic. 

In this paper we consider a two-dimensional nearest-neighbor lattice architecture using 
logical noisy SWAP gates as the basic qubit transport mechanism. It is clear that any viable 
layout for stationary qubits has to be in a two-dimensional plane so that classical control 
fields can access the qubits from the third dimension. In the case of a quasi-lD or self- 
similar architecture (such as the H-tree), the classical controls could also lie in the two- 
dimensional plane itself In our analysis, we consider the simplest possible noise model, namely 
independent adversarial stochastic noise. In this noise model, each location £ independently 
undergoes an error with some probability jg. The error that occurs can be chosen adversarially 
so that in our threshold analysis we assume that for every location the worst among the three 
Pauli errors X, Y, Z occurs. 

We find, somewhat surprisingly, that the threshold noise value for our 2D lattice architec- 
ture is close to the threshold noise value for a nonlocal architecture — both are O(10~ 5 ). We 
understand this small dependence on locality to arise from (1) the dominance of the noisy 
two-qubit gates and (2) the fact that we have optimized the circuits to be space-efficient in 
a dense, 2D square lattice, with a surprisingly small cell size of 6 x 8 for the [[7,1,3]] code, 



d A purely one-dimensional fault-tolerant architecture is impossible if one uses a distance-3 code. The reason 
is that such an architecture necessitates swapping data qubits inside a block which may generate a two-qubit 
error due to one failed SWAP. For a distance-3 code, such errors cannot be corrected. This implies that the 
level-1 error rate can never be smaller than the base error rate. A purely one-dimensional fault-tolerant 
architecture is not excluded for higher distance codes, but the threshold noise value may be very poor. 
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resulting in a small amount of moving. 

Given the good performance of our SWAP-based architecture, we believe that it is unneces- 
sary to seek alternative means of transportation (teleportation and entanglement purification 
for example), at least for the [[7,1,3]] code in this 2D layout. As the reader will observe, the 
noisy SWAP gates simply model any physical scheme in which qubits are transported through 
some quantum channel. In our model, this channel will have a fixed time-delay and a noise 
level that depends on the choice of noise level for the SWAPs. It will also be clear from our anal- 
ysis that no really long-distance means of transportation is needed for fault-tolerant quantum 
computation. 

1.1 A Local 2D Lattice Architecture 

Let us assume that we have modified an unencoded circuit M so that interactions among the 
qubits are nearest-neighbor on a 2D lattice. This requires adding SWAP gates in the unencoded 
circuit, which can be implemented using three controllcd-NOT (CNOT) gates. We will also be 
using SWAP gates in our encoded circuit, but we will treat these gates as elementary (and not 
composite) gates. 

After one level of encoding using Steane's [[7,1,3]] code, every qubit in M is replaced by a 
6x8 cell of qubits, 48 qubits in total, as depicted in Fig.Q] In a concatenated architecture, the 
next level of encoding is obtained by replacing each qubit in the cell again by a 6x8 cell, and 
so on. In this way, two-qubit CNDT gates between nearest-neighbor qubits in the cell become 
(transversal) CNOT gates between two neighboring cells. The cell functions as (1) the spatial 
unit in which error correction takes place, (2) the 'resting place' for a block of data qubits, 
and (3) the space to create a special type of ancillary encoded state \ag) (see Sections l3.2l and 
15. 3|) that is used in the computation. 

For error correction, the 48 qubits of the cell arc divided into three sets, see Fig. [TJ There 
are seven physical data qubits labeled by d in the figure. In the snapshot, the data qubits sit 
at a special location in the cell that we refer to as home base. When the data qubits interact 
with ancilla qubits for syndrome extraction, the data qubits will always be sitting at home 
base. The cell also contains at most ten ancilla qubits at a time. Three of these ancilla qubits 
are verification qubits, labeled by v in the figure. The other seven are the proper ancilla qubits 
(labeled by a) that will receive the error syndrome. The rest of the cell is filled by dummy 
qubits (labeled by O). The state of the dummy qubits is irrelevant for the computation; they 
are only used for swapping with data or ancilla qubits. One can view the dummy qubits as 
forming channels through which the other qubits can pass. Since ancilla qubits remain local 
to the cell and only data qubits interact with data qubits in neighboring cells, we place data 
qubits on the exterior of the cell and devote the interior region of the lattice cell to ancilla 
preparation and verification for error correction. 

Interactions between data qubits in neighboring cells can occur between two lattice cells 
situated on a horizontal axis, i.e., data qubits must move horizontally from home base to 
interact, or on a vertical axis, i.e., data qubits must move vertically from home base to 
interact. We include both a vertical and a horizontal transportation channel, consisting of 
dummy qubits, to be used for horizontal and vertical transportation of data qubits. The 
channels occupy the first row and first column of the lattice cell and they can be used by 
qubits in neighboring cells as well (see the discussion in Section [373)) . 
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We assume SWAP to be our mechanism for qubit transport, where a qubit moves by swap- 
ping with an adjacent qubit. To guarantee the fault tolerance of each SWAP operation we 
only swap a data or ancilla qubit with a dummy qubit. Thus, in the layout shown in Fig. [TJ 
we alternate dummy qubits with data or ancilla qubits to allow for proper fault-tolerant 
movement. 
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Fig. 1. Snapshot of the operations and state of a 6x8 lattice cell during the ancilla preparation part 
of error correction. Ancilla qubits are labeled with a and data qubits are labeled with d. Ancilla 
qubits labeled with v are used for verification. O locations represent dummy qubits that are used 
for swapping operations. The inner dashed region is the area used for ancilla preparation and 
verification. The vertical and horizontal dummy qubit border regions represent channels (labeled 
as O) used for qubit transportation. Single arrows represent CNDT gates, double arrows represent 
SWAP operations, and P z represents preparation of a |0) state. 



2 Review of Definitions and Analysis Method for the [[7,1,3]] Code 

We now review some definitions and properties of the standard fault-tolerant analysis for 
the [[7,1,3]] code; this has been described in detail in Ref. [10]. For the [[7,1,3]] code, the 
Clifford-group gates (two-qubit CNOT, singlc-qubit Hadamard (H), and single-qubit phase gate 
(S), i.e., |0) (0| + i |1) (1|) are transversal. In order to make the gate set universal, one needs 
a non-Clifford-group gate; we choose the single-qubit gate T = e ~ l7rZ / 8 . We believe that this 
choice for a non-Clifford-group gate is better than trying to implement, by local means, a 
three-qubit Toffoli gate. 

For reasons discussed later, we choose to implement the gate T, as well as that gate S, as 
in Fig. H Note that T = R z {n/4) and S = R z (n/2), where R z {9) = e~*K In these circuits, 
the state \ag) = Rz{6) |+). The corrective operation R z (29) for the T gate is the S gate and 
the corrective operation for the S gate is Z. We emphasize that with this construction, the 
only gates in the unencoded circuit are X, Y, Z, H, and CNOT (in addition to state preparations 
and measurements). 

Locations in a quantum circuit are defined to be gates, single-qubit state preparations, 
measurement steps, or memory (wait) locations; they can all be executed in a single timestep. 
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Fig. 2. Circuit for implementing the single-qubit gate Rz(0). If the measurement result is |1), then 
we correct the state by applying Rz(29). The output state is Rz(9) \ip). 



After one level of encoding, every location (denoted as O-Ga) is mapped onto a 1-rectangle 
(1-Rec), a space-time region in the encoded circuit, which consists of the encoded gate (1-Ga) 
followed by error correction (1-EC), as shown in Fig. [3] For transversal gates, the 1-Ga consists 
of performing the O-Gas on each qubit in the block(s). The encoded version of a preparation 
location is called a 1-Prep. For the [[7,1,3]] code, it consists of at most two attempts at 
preparing the desired encoded state; each attempt is followed by a verification step. If the 
verification of the first preparation fails, we repeat the preparation and verification. If it fails 
again, we continue with the failed encoded state which may contain two errors. To make a 
TRec, we follow the 1-Prep by a 1-EC. The encoded version of a measurement location is 
called a 1-Meas. In the [[7,1,3]] code, it is implemented as seven parallel measurements. In 
our model, we consider both measurement in the X basis and measurement in the Z basis to be 
single locations (rather than converting X measurement into Z measurement using a H gate). 
For the measurement 1-Rec, the 1-EC is simply the classical processing of the measurement 
outcomes. 




Fig. 3. A 1-rectangle (1-Rec), indicated by a dashed box, that replaces a single-qubit O-Ga location. 
The 1-Rec consists of the encoded fault-tolerant implementation of the O-Ga (1-Ga) followed by 
an error-correction procedure (1-EC). 



For the fault-tolerance analysis, one also defines an extended 1-Rcc (1-exRcc) which con- 
sists of a 1-Rec with its preceding l-EC(s) on the input block(s). The 1-exRec for a preparation 
or measurement location coincides with the 1-Rec for preparation or measurement. 

In order to be fault-tolerant, the 1-EC, 1-Ga, and 1-Prep must satisfy the fault-tolerance 
requirements listed in Ref. jTU] and reproduced in Appendix [TJ Our constructions obey these 
requirements with the caveat that the 1-EC does not satisfy the second property. As we will 
show, it can leave an X error and a Z error on the data after its application. Thus for fault 
tolerance, it is necessary that the 1-Gas never transform a X or Z error into a Y error since this 
would correspond to having two X (or Z) errors in a block. This is achieved by only having 
O-Gas which are X, Y, Z, H, and CNOT (and not S). 
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If the fault-tolerance requirements are fulfilled, it can be shown that if a 1-exRcc is good 
— for a distancc-3 code, such as the [[7,1,3]] code, this means that the 1-exRec contains at 
most one fault — then the dynamics induced by the 1-Rec contained in the 1-exRec is correct 
[10] . as defined in this picture: 



correct 
1-Rec 



perfect 
1-decoder 



perfect 
1-decoder 



ideal 
O-Ga 



A 1-cxRec containing two faults can be bad, but not all higher-order faults lead to two 
errors on the data. In order to give an accurate estimate of the threshold, Ref. [10] introduced 
the notion of benign and malignant fault patterns. A set of locations is benign in a 1-exRcc 
if the 1-Rec contained in the 1-exRcc is correct. Otherwise, the set of locations is called 
malignant. With this definition one can redefine a bad or failed 1-exRec, namely, as a 1- 
cxRec that contains faults at malignant locations. In our analysis, we assume any pattern of 
three faults or more is always malignant. We discuss this assumption in Section 14. II since it 
has a small but non-zero effect on the threshold estimate. 

We denote the failure probability of a location of type £ as je = 7°. The failure probability 
at encoding level r for location £ can be computed with a multi-dimensional map F, with 
Ft (7 r ~ 1 ) = 7J. The map F is independent of concatenation level and is determined by the 
number of malignant and higher-order errors in the 1-exRec for the different locations. An 
initial vector of probabilities 7 is said to be below the threshold if there is a concatenation 
level r e such that 

Ve,3r e ,W, 7? = (F^( 7 °)) £ < e. (1) 

In Sections |4] and El we describe how we calculate this map and the threshold results that 
we obtain for the 2D local and the nonlocal model. We now describe the local fault-tolerant 
circuits that we use in our analysis. 

3 Local Fault-tolerant Circuitry 

Before we describe our error-correction circuitry, we first describe local 1-exRecs. In a local 
architecture, a 1-Ga which replaces a two-qubit 0-Ga includes the necessary SWAP operations 
to move the two encoded qubits together. At the next level of concatenation, each new 
location, including a SWAP, is replaced by its 1-Rec. Our SWAP gates involve one dummy qubit 
and one data qubit; it is unnecessary to do error correction on the dummy qubit. Therefore, 
the SWAP 1-Rec effectively behaves as a 1-Rec for a single-qubit gate since it requires a 1-EC 
on only one block of qubits. 

3.1 Local Steane Error Correction (1-EC) 

For [[7, 1,3]] error correction, we have a choice between (1) Shor-type error correction, where 
the syndrome for each stabilizer operator is measured separately, (2) Stcanc-type error cor- 
rection, where the entire X- or Z-error syndrome is copied onto a verified ancilla block, or 
(3) Knill-type error correction, where the data qubits arc telcported into a new block and 
X- and Z-error-syndrome information is obtained via a Bell measurement. For our local lay- 
out, it is important to minimize the use of space and thus movement, potentially at the cost 
of increased waiting and memory errors. The reason is that we assume memory errors will 
typically occur at lower rates than SWAP errors. This optimization would seem to point to a 
sequential execution of the most space-efficient error correction, which is Shor's. However, in 
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Shor-type error correction, the number of gates between data and ancilla qubits is more than 
in Steane-type error correction, which can negatively affect the threshold. 

Therefore, we opt for a variant of Steane's error correction, schematically drawn in Fig. [4] 
We use a deterministic error-correction network, in which we sequentially prepare three ancilla 
blocks in the encoded 0-state |0) . If a single failure causes one ancilla block to fail verification, 
then two ancilla blocks remain, one each for X- and Z-error correction. If two ancilla blocks 
fail verification, which by construction can only happen when two faults occur, we cannot 
perform full error correction and we count that as a failure of the data in this rectangle. If 
the first two ancilla blocks pass verification, then there is no need to prepare a third ancilla 
block. Instead, the data block waits for the same number of timesteps that would have been 
required for the third ancilla block preparation. 
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Fig. 4. The error-correction network (1-EC) based on Steane's [[7, 1, 3]] code, performing first 
Z-error correction, and then X-error correction. Error correction consists of ancilla preparation 
((/), ancilla verification (V), syndrome extraction (5), and a recovery operation (TV). In many 
circumstances, the application of the corrective gate 71 can be postponed and is carried forward 
classically as a change of 'Pauli frame' [7|. 

The local versions of Q and V are shown in Fig. [5] Our ancilla verification scheme is less 
extensive than the standard CSS method described, for example, in Ref. [10]. In the standard 
method, the ancilla is verified using another block in the state 1 0) . For the fault tolerance of our 
error-correction network, we only need to measure the single stabilizer IZIZIZI on the ancilla 
block |0) f The ancilla needs to be verified since the encoding uses two-qubit gates which 
can put two errors on the block. If these errors are X errors they can propagate to the data. 
However, not every possible pair of X errors will be generated by the encoder circuit. In fact, 
we only need to check for errors on the second, fourth, and sixth qubits of the ancilla block. 
Our verification method prevents an ancilla from passing with the errors X2X7, X3X6, X4X5, 
each of which can be caused by a single faulty location in the preparation circuit. An ancilla 
with single X and/or single Z errors can still pass verification. To parallelize the necessary 
measurement, verification is performed by preparing a 3-qubit cat state (|000) + | lll))/v / 2 
and applying controlled-Z (CZ) gates between the three verification qubits and the second, 
fourth, and sixth qubits of the ancilla block. Each verification qubit is then measured in the 
X basis. If the parity of the measurement results is 1, the ancilla block fails verification. If 



e Thanks to Panos Aliferis for suggesting this simplification. 
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an ancilla block passes verification, it is used for syndrome extraction (5), shown in Figure [6] 
During ancilla preparation the data qubits are moved to home base in order to interact with 
the prepared ancilla. 
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Fig. 5. The local preparation and verification networks Q and V. The lower three qubits are 
the verification qubits of network V. The seven timestcps in which this circuit is executed are 
delineated by the long vertical lines. The third timestep is shown in Fig. [l] The single-qubit X 
symbol indicates that in this timestep, the qubit has been swapped with a dummy qubit. The CNOT 
symbol with the control indicated with a X means a CNOT followed by SWAP on the same qubits. 
The two-qubit gate indicated by a solid circle on one end and a X on the other is a CZ followed by 
SWAP. It is important to note that the qubit labels on the left correspond to the associated wires 
throughout the circuit (and thus are never swapped) and thus qubit a6, for example, interacts 
with qubit v3 in the 5th timestep. Thus the fourth qubit travels the farthest, being swapped four 
spaces along the square lattice during the operation of this circuit. 



3.2 Local Preparations (1 -Preps) 

The local 1-Prep for |0) consists of at most two attempts to create a verified ancilla as in the 
local network in Fig. O We obtain the 1-Prep for |+) by applying transversal Hadamards on 
the passed ancilla. For universality, we also need 1-Prcps for the state \a$), for = 7r/2,7r/4 
(see Fig. [5]). We can produce such states as follows (see AGP [TO]: The unencoded state 0^/2) 



is correctly generated by measuring the operator SXS^ on the state |0), since SXS^ 



1 ir/2 



a 7r/2)- If wc nn d eigenvalue +1, we arc done, otherwise we apply Z. Similarly, for 0^/4), 
we measure TXT/t and correct with Z if eigenvalue —1 is found. As we will see later, \ag) 
preparations are fundamentally different from Clifford-group preparations, in that it is not 
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Fig. 6. The network S for syndrome extraction during (a) X-error correction and (b) Z-error 
correction. During syndrome extraction the data qubits sit at home base in the cell. The two- 
qubit gate with solid circles on both ends is the CZ gate. 



necessary to make them fault tolerant at every level of concatenation. 
3. 3 Local Rectangles and Their Use of Space and Time 

Table 1. Types and number of locations contained in a local 1-EC and two local 1-Gas. The 
first numbers for the 1-EC are based on the assumption that the first and second ancillas pass 
verification; the numbers in parentheses result from the first or second ancilla failing verification. 
The last column lists the number of locations in a horizontal CNOT 1-Ga. 
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In Table[TJ we list the universal set of location types considered in our scheme. Preparation 
of |+), measurement in the X basis, the CZ gate, and CNOT and CZ combined with SWAP are 
considered single locations which can be implemented in one fundamental timestep. We will 
assume that these composite gates have error probabilities similar to the other elementary 
gates. In practical implementations, some gates may have to be constructed from multiple 
gates and may thus have higher error rates. We assume that the time it takes to do a 
measurement is the same as the single-qubit gate time and that classical post-processing docs 
not take any additional time. Note in Tablc[T]that we use two different types of wait locations, 
one that happens in parallel with a gate and one that happens in parallel with a measurement. 
The reason for having separate locations is that the encoded version of the wait location that 
occurs in parallel with a gate needs many additional wait locations in order to be synchronized 
with the encoded version of a gate. 

We include both horizontal and vertical two-qubit gate locations. Two-qubit gates occur- 
ring between two logical qubits x and y on a horizontal axis use the horizontal channel internal 
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to lattice cell x or lattice cell y as well as the horizontal channel internal to the lattice cell 
belonging to the neighboring qubit directly south of lattice cell x or lattice cell y (i.e., directly 
below it). We will actually use either the channel south of lattice cell x or the channel south 
of lattice cell y, but not both. For example, in Fig. Q^a), data qubits d6, d5, and d3 move 
into the channel internal to their lattice cell, while qubits dA, d2, dl, and dl move into the 
neighboring channel directly south. To perform a two-qubit gate, the neighboring horizontal 
qubits then move across the lattice cell by swapping (for example, along row 2 and 6 in the 
figure) while the qubits that moved into the channels remain stationary and wait. Vertical 
two-qubit gates only use the vertical channel internal to the lattice cells of the two qubits. 
It follows that two vertical two-qubit gates on two adjacent pairs of cells can be executed 
in parallel since there is no common use of channels. For two horizontal two-qubit gates on 
pairs of cells that are above each other we need some scheduling since a common channel 
is used. This scheduling is easy to achieve, since we can designate which of the two qubits 
involved in the two-qubit gates will move into the channel, and which will remain within the 
cell. The qubits remaining within the cell can swap across the internal row freely, while the 
qubits moving into the channel remain stationary and wait until the data qubits moving along 
the interior of the lattice cell arrive. For a more detailed explanation, a complete snapshot of 
the movement just described is available at 

http://www.research.ibm.com/quantuminfo/svore/movies.html . 

In the local model, we must clearly define to what rectangle some of the SWAP operations 
and additional wait locations belong since this can be ambiguous. During ancilla preparation 
in the 1-EC, the data qubits move around in the cell and reside at home base once syndrome 
extraction begins. This movement and possibly waiting of data qubits is included in the 
definition of the previous 1-Ga, since it can be viewed as movement back to home base after 
the 1-Ga was performed. Similarly, the 1-EC does not contain any SWAP or wait locations on 
the data once the last syndrome extraction (and possibly error correction) has been performed. 
Thus the 1-Ga includes all data locations between the last time the data was in home base 
for syndrome extraction and the next time the data is in home base for syndrome extraction. 
The 1-EC includes all locations on the data between the first and last syndrome extraction, 
including these syndrome extractions (and possibly error correction) themselves. 

With these definitions, we list the number of locations in a few essential routines in Table 
[U Note that the total number of SWAP gates in the h-CNOT 1-Ga is 112 + 14 = 126 which is 
substantial. As we will see in Section [5] the effect on the threshold of these additional gates 
is quite small. 

In Table [2j we list the number of timesteps and locations in a 1-Rec and 1-exRec for 
different locations £. The local 1-EC takes 27 timesteps; this may be compared to the same 
1-EC without SWAP, i.e., the nonlocal 1-EC, which takes 21 timesteps. Note that in Table 
[2] the difference between the number of timesteps in a 1-Rec and a 1-exRec is 26, which is 
one less than the number of timesteps in a 1-EC. This is because there is an overlap of one 
timestep between the leading 1-EC and the 1-Rec. 

Of all gates the horizontal CNOT 1-Rec takes the longest, a total of 35 timesteps (compared 
to 22 timesteps in the similar nonlocal model). The single-qubit 1-Rcc takes 28 timesteps 
(compared to 22 timesteps in the nonlocal model). For synchronization, we have added wait 
locations to the 1-Rccs that take fewer than 35 timesteps, so that the time for all 'gate' 1-Recs 
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is 35 timcstcps. Notice that the preparation 1-Recs take more time, but these 1-Recs can be 
started at the appropriate time in advance so that the prepared states will be ready when 
necessary. 

Since it is hard to represent what happens in three-dimensional space-time on a two- 
dimensional sheet of paper, we only partially represent in this paper the local fault-tolerant 
circuits that we have developed. Supplementary material in the form of complete sequences 
of time snapshots of the 2D layouts can be found at 

http://www.research.ibm.com/quantuminfo/svore/movics.html. 

Tabic 2. The number of timesteps in the 1-Rec and the 1-exRec for different locations £. 



e 


Time of 1-Rec 


Time of 1-exRec 


Total # loc. in 1-exRec 


[0-7] 


35 


61 


1225 


[8-10] 


35 


61 


616 


[11-12] 


41 


41 


469 


[13-14] 


1 


27 


196 


[15] 


35 


61 


616 


[16] 


28 


54 


378 



4 Methodology 

Our threshold analysis relies on counting the number of malignant fault pairs in a 1-exRec; 
thus it is important to understand this notion. In principle, we have to assume an arbitrary 
input to the 1-exRec in order to analyze malignancy of faults. Given the fact that the [[7,1,3]] 
code is a perfect code, incoming errors can be modeled by letting the input block have at most 
1 X and 1 Z error. Now one considers two Pauli faults in the 1-exRec. If both faults occur in 
a single leading 1-EC, the output of that 1-EC can again be viewed as some codeword plus 
at most one X and Z error. This implies that the subsequent 1-Rec is correct, and thus the 
fault-pattern is benign. If the 1-cxRec is a CNOT 1-exRcc, both leading 1-ECs can have faults 
and the transversal CNDT gate could spread these so that the 1-Rec may not be correct. 

Let us consider how the possible input errors affect the output of a 1-EC with at most 
one fault. If the 1-EC has no fault, the incoming errors will be corrected and the output 
is a codeword, just as if there were no incoming errors. If the 1-EC has a single fault, it 
is guaranteed that one of the incoming errors (X or Z) will be corrected. As was argued in 
Rcf. [10] , the overall effect of the other incoming error and the fault in the 1-EC is a possible 
logical operation in the code space plus an error whose position and character only depends 
on the fault in the 1-EC. In other words, the errors at the end of the leading 1-EC with at 
most one fault do not depend on the pattern of incoming errors. This property still holds 
with our modified 1-EC where the 1-EC can create one Z and one X error in the outgoing data 
block. This implies that for determining malignancy, we only need to place two faults inside 
the 1-exRec and not consider incoming errors. Since we are using an adversarial noise model, 
we consider a pair of locations malignant if it is malignant for some Pauli errors occurring at 
those locations. 

4-1 Calculation of Failure Probabilities 

We calculate the failure probability 7J = F^(7 r_1 ) of a type-£ 1-exRcc. Note that F is 
level-independent and we drop the superscript r in what follows. We can bound the failure 
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probability 7^ as 

P(failurc of a 1-cxRcc) = P(l-Rcc is incorrect) 

< P(l + trailing 1-ECs do not occur) 

+ P(l-Rcc is incorrect & trailing 1-ECs occur). (2) 

Remember that a 1-EC does not occur when less than two ancilla blocks pass verification. 
Consider P(l-Rcc is incorrect & trailing 1-ECs occur). Because of the fault tolerance of the 
quantum circuits, the 1-Rcc is only incorrect when at least two malignant faults occur in the 
1-exRec. We can upper-bound this probability as 

P(l-Rec is incorrect & trailing 1-ECs occur) < P(l-Rcc incor. due to mal. pair 

& trailing 1-ECs occur) 
+ P(3 + faults & trailing 1-ECs occur 1(3) 

Consider the failure probability due to malignant pairs. In this probability it is not specified 
whether the leading 1-ECs in the 1-exRec occurred or not. But the only way for a leading 
1-EC to fail is when there are two faults in it. As argued earlier in this section, for the [[7,1,3]] 
code it follows that the 1-Rec is correct and thus we can safely assume that the leading 1-ECs 
always occur. We estimate this probability using a numerical malignancy count 

P(l-Rcc incor. due to mal. pair in 1-exRcc of type £ 

& all 1-ECs occur in 1-exRcc) < V^ali uffilj- (4) 

i>j 

Here of, r 1 is the number of malignant pairs of type i and j for a 1-exRcc of type £, where 
we count only those malignant pairs that do not lead to failure of a 1-EC. The numerical 
malignant pair counts are obtained using the QASM Toolsuite [23). With this software, one 
simulates the propagation of Pauli errors in a quantum circuit and thus determines whether 
a pair of faults is malignant. A selection of these malignancy matrices afi r ■-. are reproduced 
in Appendix [Bj 

We upper-bound the triple-fault term as P(3 + faults & trailing 1-ECs occur) < (3)7^, 
where N is the total number of locations in the full 1-exRec (including all 1-ECs) and 7 max = 
max^(7^) is the maximum failure probability. This expression is probably a large over-estimate 
of the higher-order terms. For our largest local 1-exRec, the CNDT 1-exRec, N = 1225. 
This implies that, including only these triple malignancies, the threshold cannot be above 
p = 5.7 x 1CP 5 (which is the solution of the equation p = p 3 ( 3 )). Since the third term is a 
very steep function of p, the effect of the term around the thresholds that we find below is not 
so large. In the results in Section[5l we see that the effect of omitting this term altogether gives 
similar thresholds in the O(10~ 5 ) range. Alternative methods of bounding this higher-order 
term exist, but they either involve triple malignancy counting (which is very time consuming) 
or heuristic strategies such as malignancy sampling. 

Let us return to the first term in Eq ((2]). We calculate the probability that one or more 
of the 1-ECs fail as 

P(l+trailing 1-ECs fail) = 1 - P(1-EC occurs) 9 , (5) 
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Table 3. Number of bad locations of type I in the Q and V networks that make the ancilla not 
pass verification. Location types not listed are all good. 



Description of type t 


Number 


h-CNOT 


1 


v-CNOT 


7 


h-CNOT followed by SWAP 


2 


h-CZ followed by SWAP 


2 


h-SWAP 


3 


|0)-preparation 


2 


+)-prcparation 


1 


X-mcasurcmcnt 


3 


wait during gate/prep 


2 



where g = for a measurement 1-cxRcc, q = 1 for singlc-qubit 1-cxRecs (including SWAP and 
1-Prcp), and q = 2 for the two-qubit 1-exRccs. Furthermore, 

P(1-EC occurs) = P(anc. pass) 2 + 2P(anc. pass) 2 P(anc. fail). (6) 

How do we estimate the probability for an ancilla to pass? We have P(anc. pass) = 
1 — P(anc. does not pass), where we bound 

P(anc. does not pass) < P(no pass due to 1 fault in Q,V) + P(2 + faulty loc. in Q 1 V), (7) 

where P(2 + faulty loc. in Q, V) < ( 2 9 )7max- Thus we assume any two faulty locations in 
preparation and verification cause the ancilla to fail verification. This is a slight overcount, 
but the majority of single faulty locations in fact cause the ancilla to fail verification. We 

have 1 



P(no pass due to 1 fault in Q, V) = \^ otfyt, 

£eG,v 



(8) 



where ae is the number of bad locations of type £ occurring in the Q, V network as determined 
by the QASM toolsuite, listed in Table [3 

5 Threshold Studies 

We determine a lower bound on the noise threshold for each location type using the formulas 
for failure probabilities given in the previous sections. In order to estimate thresholds, we 
repeatedly apply the 17-dimensional flow map ^(7) to the initial vector of failure probabili- 
ties. The only flow that depends on the a„/2~ and a ff /4-preparations is the flow map for these 
preparations themselves. This implies that the threshold for all the other (Clifford) gates is 
independent of the initial failure probabilities for these two preparations. Thus we can esti- 
mate the threshold for these preparations separately and see whether or not they determine 
the overall threshold. 

5.1 Clifford-group Location Thresholds 

Let us first consider a much studied case in the literature in which all initial probabilities are 
the same, i.e., set to some value 7, except the failure probability for a memory location which 
is set to 7/IO. The threshold of all these gates (except the memory location) will be denoted 
as 7 C . For our 2D lattice architecture, we find a lower bound of 7 C > 1.85 x 10~ 5 . We compare 
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our local 2D architecture with a nonlocal architecture with exactly the same circuits (except 
for the SWAPs). For the nonlocal setting, we find j c > 3.61 x 10~ 5 . 

In Section 14.11 we argued that we considerably overestimate the three-or-more- fault terms 
in our failure probability calculation. To see how this affects the threshold, we perform our 
analysis omitting this triple- fault term entirely in Eq ([3]). The local threshold then becomes 
2.15 x 10~ 5 and the nonlocal threshold becomes 4.32 x 10~ 5 , which are small corrections on 
the numbers given above. 

To test the correctness of our local analysis, we also analyze our local 2D architecture in 
a variant where at all levels of concatenation the error rates of SWAP gates are set to zero by 
hand. Then the only difference between the local and nonlocal case is the additional wait 
locations in the local architecture. In this scenario, the threshold is 7 C > 3.0 x 10~ 5 , in which 
we see the effect of additional waiting in the local model. To eliminate the possibility that 
the threshold for the nonlocal architecture is suboptimal due to the choice of deterministic 
error correction, we also analyze the nonlocal circuits in [10J : this gives a lower bound of 
7c > 4.3 x 10~ 5 . This shows that our choice of circuits is close to optimal. 

One may think that the closeness of the local and nonlocal thresholds could be due to the 
fact that we let the memory error rates be 1/10 of the gate error rates and thus eliminate 
the effect of additional waiting on the local threshold. Thus we also consider the local and 
nonlocal architecture where all error probabilities are identical. For the local architecture we 
obtain a lower threshold, -f c > 1.1 x 10 -5 ; for the nonlocal architecture the threshold becomes 
1.9 x 10~ 5 . 

To understand the behavior of the failure probabilities at different levels of encoding we use 
a threshold reliability information plot (TRIP) [9]. The behavior of the failure probability for 
the horizontal CNOT in the local model is given in Fig. [7l The value of 7 at the crossing point 
with the 45-degree line is the threshold at level r, also called the level r pseudothreshold, for the 
given location type. In Figure[51 we have plotted the behavior of the SWAP location. From these 
figures it is clear that the horizontal CNOT dominates the threshold. The pseudothrcsholds of 
the SWAP location get worse as a function of r due to the noisiness of the CNOT, whereas the 
pseudothresholds of the CNOT get better as a function of r. Another way of seeing this is that 
after one level of encoding there is a region in which one is below the pseudothreshold of the 
SWAP gate, but not below the pseudothreshold of the CNOT gate. In this region, the CNOT gets 
worse, which feeds into a worse pseudothreshold of the SWAP at the second level of encoding. 

5.2 Approximate Closed Form for Thresholds 

Even though our map is high-dimensional, we can capture the essential features by a two- 
dimensional map. We do this by making some drastic approximations: the 1-EC failure 
probabilities of Eq. ((2|), and the 3+ fault probabilities of Eq. (J3j> , are ignored. Many distinct 
rectangles are treated as having the same error probabilities. Despite these approximations, 
the thresholds can be estimated very simply to about a factor of two. 

We assume that all error levels of gates, preparations, and measurements are comparable 
(but not necessarily identical) at the base level. If some of the error rates are much larger 
than others (say, by a factor of 100), it may take more levels of concatenation for the two- 
dimensional representation to become a good approximation, but these initial conditions will 
be eventually marginalized. 
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Fig. 7. Initial failure probability 7 versus the failure probability 7J at level r of a horizontal CNOT 
location, shown for levels r = 1, . . . , 4. 
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Fig. 8. Initial failure probability 7 versus the failure probability at level r of a horizontal SWAP 
location type 7J , shown for levels r = 1, . . . , 4. 



We only consider the Clifford gates. After a single iteration, information about the loca- 
tions in each 1-exRec is absorbed into the new error probabilities of simulated locations. Then 
we find if we start with comparable error rates for all gates, all subsequent iterations have the 
property that the failure probabilities fall into three distinct groups: (1) those of two-qubit 
gates (locations [0-7]) (2) those of single-qubit gates (locations [8-10], [15], [16]), (3) those 
of the preparations (locations [11], [12]) and those of measurement (locations [13], [14]). In 
other words, the failure probabilities of members of each group are very close in value. We 
see that this grouping basically corresponds to the number of 1-ECs in the 1-cxRcc. 

We also observe that the measurement failure probabilities are very small compared to 
the other failure probabilities and that the behavior of the map changes very little if the 
measurement failure probability is always set to zero. After the second iteration, the failure 
probability of the preparation rectangles gets close to those of the single-qubit gates, so 
we group them together. Thus, we end up with a two-dimensional map described by the 
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equations: 



x i+1 



»'« = (x' „<) • ( f ;)•(•). (io) 

where i is the iteration level. Here x l+1 is the failure probability of the two-qubit gate group 
at level i + 1 and y l+1 is the failure probability of the singlc-qubit gate group. The numbers a, 
b, etc. can be extracted by summing entries in the malignancy matrices af , r -, (see Appendix 
[B]for two examples of single-qubit and two-qubit gates). For the local model, the numbers for 
£ = (horizontal CNOT) and I = 9 (horizontal SWAP) are a = 7907, b = ^1 and c = 93488 
and d = 1956, e = ^|M and / = 35886. 

The equations for the fixed point of this map are two coupled second-order equations 
for x and y. By elementary elimination methods, one can derive from these a fourth-order 
polynomial in x whose roots give the fixed point. One root is at x — 0; factoring this out 
leaves a cubic polynomial. We find that, for the range of parameters we have, the next root, 
giving the threshold for x, is accurately given by truncating this cubic polynomial to linear 
order. The resulting expression for the threshold value of x is 

(11) 



ac + Ace - 2b f - ,/ 2 ' 
A parallel analysis for y gives an expression for its fixed point: 

V = a 2 + Abd - 2ae + df ' ^ 2 ' 

Now we analyze the first iteration of the full map. For simplicity, we try to follow the 
analysis in the previous section where all error rates are the same, except for the memory 
error rates. At the physical level, we start with a gate failure rate 7 fl and a wait failure rate 
j w . The a; 1 and y 1 are given by quadratic expressions 



.?■ 



'■w-iJ2ili' (13) 



»'-<■**>■(**)■(£)■ (14) 

For the local model, the parameters arc F = 71779, G = ^|^, H = 26842 (taken from the 
horizontal CNOT malignancy matrix) and J = 8318, K = ^|43 anc j £ _ 21632 (extracted from 
the singlc-qubit wait malignancy matrix). 

We have been interested in the case *y w = ^7. We find this case is numerically similar 
to setting j w = 0, causing these equations to become very simple: x 1 = -F7 2 , and y 1 = J7 2 . 
Setting x 1 = x, the threshold expression, and y 1 = y (the threshold expressions Eqs (|11I12[) 
and solving for 7, we get two different expressions for 7: 



7 x > F{ac + Ace-2bf- p) 1 ^ 
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7 y J(a 2 + Abd - 2ae + df) ' ^ 

We will be below threshold if 7 is lower than either of these two expressions, since then both 
x 1 and y 1 will be below their threshold values. Therefore, 



7c mm \^J f(qc + ^ _ ^j _ py \j J(fl2 + m _ 2ag + d f) J ' ^ 

Note that these expressions satisfy the required homogeneity property, namely that if all 
malignant pair counts are scaled by the same constant k, the threshold changes by the factor 
1/k. Putting in values for these malignant pair counts, we find that the estimated threshold 
values are 7 C > 1.94 x 10 -5 for the local model and j c > 4.67 x 10~ 5 for the nonlocal model. 
This agrees very well with the threshold values of the actual approximate map, which arc 
1.94 x 10~ 5 and 4.43 x 10 -5 ; these are in turn in reasonable agreement with the threshold 
values for the exact map, which were found to be 1.85 x 10 -5 and 3.61 x 10 -5 . 

Our formulas explain the factor-of-two difference between the local and nonlocal cases in 
the following way. The malignant pair counts differ between the two models, in different cases 
by a factor of 5 at most, and a factor of 1 (unchanged, that is) at the least. These formulas 
represent a kind of weighted average among these different counts, so the actual difference is 
somewhere between a factor of 5 and a factor of 1 . 

5.3 ag -preparation Threshold and Injection by Teleportation 

It would be possible to approach ag-preparation in exactly the same way as all other locations 
in the quantum circuit. A deterministic, fault-tolerant extended rectangle for this preparation 
is known, which is based on the fact that the state \ag) is an eigenstate of the encoded 
operator R z (0)XR z (—0) = R Z (28)X. Fault tolerance is achieved by repeated measurement of 
this operator followed by error correction, exactly as in Fig. 13 of AGP [10]. Since we require 
this rectangle to be deterministic, we must allow for three repetitions, rather than the two 
repetitions used in AGP. 

Unfortunately, this triple-repetition cxRcc, with all the moves and waits required by local- 
ity, has an extremely large number of locations, and the threshold, if we determine it using the 
formalism discussed for the Clifford gates, is very poor. But we can take advantage of the fact 
that the \ag) preparations are not mandatory at the lowest level (because these preparations 
are not needed in 1-ECs) to obtain fault tolerance by a different strategy. 

This strategy is injection by teleportation, an idea due to Knill [7j. The key circuit is 
shown in Fig. [9l here, an encoded EPR pair is created and one half of the pair is decoded. 
This decoded half is used to teleport an unencoded \ag). After teleportation, the resulting 
state is the encoded \ag). The encoded state produced by this circuit may be at any level of 
concatenation. This circuit is not fault-tolerant, so \ag) may be noisier than \ag). But this 
circuit will only suffer from single faults in gates in the latter part of the decoding circuit 
and in the measurements. The total number of such gates is on the order of ten, no matter 
what the level of concatenation of the final, coded state. Thus, if we inject to level k of 
teleportation, the encoded error rate 7^ is given approximately by 7^ = 7^ + 107° NOT . 

For purposes of analysis, wc imagine that the injection level k is chosen so that the Clifford- 
group gates are essentially noiseless. At this level, we may use a fault tolerant circuit related 
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to Fig. 13 of AGP [10], as described above. Upon further concatenation, the probability 
flows will be determined just by the preparation counts in this 1-exRcc. There would be 84 
a,0-prcparations in the triple-repetition rectangle; the ag-thrcshold, then, is estimated by the 
equation jg = ( 2 )le- Assuming 7 CNO t = 10~ 5 , these equations give an estimate for the jg 
threshold of 1.0 x 10~ 4 . 

Thus, the a^-threshold value is not unduly low; we expect, in fact, that this number can be 
improved with modifications of the AGP circuit. We note that this threshold is of course far 
from the theoretical upper maximum of about 17% established by the analysis of the magic- 
state distillation technique [24). We also note that, with the scheme just outlined, there is 
a serious efficiency cost, in that more levels of concatenation are typically needed to achieve 
effectively noiseless operation. We expect that this, too, may be improved upon. 
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Fig. 9. Injection by teleportation. An encoded EPR state is created and half of the state is 
decoded. A Bell measurement is carried out on the state \ag) and the decoded half followed by 
encoded corrective Paulis on the encoded half of the EPR pair. There is ample room in the two- 
dimensional cell to lay out this circuit, whatever the level of concatenation desired for the final 
state. 



6 Concluding Remarks on the Noise Threshold of the [[7,1,3]] Code 

In the fault-tolerance literature, various estimates of the nonlocal threshold for the [[7,1,3]] 
code have been given. These numbers vary between an optimistic value of O(10 -3 ) [5], an 
estimate of O(10 -4 ) [5] and the most recent rigorous lower bound of 2.73 x 10 -5 established 
in Ref. [lOj . In this section, we argue why, for stochastic independent noise models, the noise 
threshold of the [[7,1,3]] code is realistically O(10~ 5 ). 

The main reasons that the earlier estimates were more optimistic lies in the method of 
analysis. In Ref. |10| one estimates the failure probability of an extended rectangle, which 
includes error correction before and after the encoded gate. The previous work analyzed 
the failure probability of a rectangle and tried to take into account incoming errors into the 
rectangle in a heuristic fashion. The reason that the analysis in Ref. [TU] is superior to the 
other analyses is that the definition of failure introduced in Ref. [10| has a direct interpretation 
in terms of the failures in the unencoded quantum circuit that the noisy quantum circuit 
simulates. Namely, by pushing perfect decoders from the end of the computation through to 
the front, one generates the unencoded circuit in which failed 1-exRec corresponds to failed 
elementary gates. 

In Ref. [TU], the malignancy method was used to establish a lower bound of 2.73 x 10 -5 for 
the [[7, 1, 3]] code in the presence of independent adversarial stochastic noise. This derivation 
of the lower bound assumed, for simplicity, that all gates in the circuit are (worst-case) 
CNOT gates. One may thus think that such an analysis would be far too pessimistic. This 
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conclusion is falsified by two observations. The first observation is that there are already 13245 
malignant CNOT pairs alone in the CNOT f-exRec (see the malignancy matrix in Rcf. [10]). 
Thus just assuming that all gates arc noiseless except the CNOT would set a threshold of 
13 J, 45 w 7.55 x 10 -5 . The second observation is that the full nonlocal analysis of the present 
paper that takes into account all types of locations similarly produces a O(10~ 5 ) threshold. 
The threshold also does not change very much if we change adversarial noise to depolarizing 
noise as shown in Ref. [10]; their lower bound for depolarizing noise is again O(10~ 5 ). 

Another cause for potential looseness in threshold estimates using the malignancy method 
is how one deals with incoming errors in the 1-exRec. For general codes, the patterns of 
incoming errors into the 1-exRec may determine which pattern of faults is malignant inside 
the 1-exRec. If we assume a worst-case pattern of incoming errors (instead of some steady- 
state pattern), our threshold estimate could be lower than necessary. But for the perfect 
[[7, 1, 3]] code, one can argue (as was done in Ref. [10] and repeated here) that the malignancy 
of faults inside the 1-exRec is independent of the pattern of incoming errors. 

These observations lead to two conclusions that are of interest in further studies. Firstly, 
it is highly desirable to look at different code architectures and see whether one can obtain 
a threshold in the range O(10 -3 ) — O(10 -4 ). One example of such a promising architecture 
is the C^/Cq scheme by Knill [7]- Secondly, as witnessed by the analysis of the [[7,1,3]] code, 
thresholds have certain crucial dependencies but are otherwise remarkably robust against 
small variations (in the encoding circuitry, in the noise rates of some gates, and in adding 
SWAP gates for locality, for example). To identify these crucial dependencies, by studying 
a whole range of codes and understanding in what part of 'code-space' the most promising 
codes lie, will be an important task ahead. 
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Appendix A Fault-tolerance Requirements 

The [[7, 1, 3]] code is a perfect code which means that all states in the 2 7 dimensional space 
can be represented as some codeword with at most 1 X error and at most 1 Z error. 

The fault-tolerance properties that the constructions of 1-Gas, 1-Meas's, 1-Prcps, and 1- 
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Table A.l. The six stabilizers and logical operations of the [[7,1,3]] code. 



Name 


Operator 


si 


IIIXXXX 


S2 


IXXIIXX 


S3 


XIXIXIX 


s 4 


IIIZZZZ 


S5 


IZZIIZZ 


S6 


ZIZIZIZ 


X 

z 


XXXXXXX 

zzzzzzz 



ECs must satisfy for a distance-3 code, as stated in Ref. [TO], are: 

1. If a 1-EC contains no fault, it takes any input to an output in the code space. 

2. If a 1-EC contains one fault, it takes any input to a valid output. (The output of a level-1 
block is 'valid' if it deviates from the code space by the action of a weight- 1 operator). 

3. If a 1-EC contains no fault, it takes an input with at most one error to an output with no 
errors. 

4. If a 1-EC contains at most one fault, it takes an input with no errors to an output with at 
most one error. 

5. If a 1-Ga contains no fault, it takes an input with at most one error to an output with at 
most one error in each output block. 

6. If a 1-Ga contains at most one fault, it takes an input with no errors to an output with at 
most one error in each output block. 

7. A 1-Meas with no faults applied to an input with one error agrees with an ideal measure- 
ment. 

8. A 1-Mcas with at most one fault applied to an input with no errors agrees with an ideal 
measurement. 

9. A 1-Prep with at most one fault produces an output with at most one error. 



Appendix B Selected Malignant Pair Matrices 



In these matrices the columns and rows corresponding to location i are indicated by [i] 
Columns (and rows) with only zero elements are omitted. 



For the horizontal SWAP, location 9, wc have 
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( [0] [1] [3] 

16 186 68 

467 312 

63 
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For the horizontal CNOT, location 0, we have 
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a [i\,\j] 



[0] 
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